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I. INTRODUCTION 


Research Institute for Advanced Computer Science 

JOSEPH OLIGER 
Director 

The Research Institute for Advanced Computer Science (RIACS) was established by the Universities Space 
Research Association (USRA) at the NASA Ames Research Center (ARC) on June 6, 1983. RIACS is pri- 
vately operated by USRA, a consortium of universities with research programs in the aerospace sciences, 
under a cooperative agreement with NASA. 

The primary mission of RIACS is to provide research and expertise in computer science and scientific com- 
puting to support the scientific missions of NASA ARC. The research carried out at RIACS must change its 
emphasis from year to year in response to NASA ARC’S changing needs and technological opportunities. A 
flexible scientific staff is provided through a university faculty visitor program, a post doctoral program, 
and a student visitor program. Not only does this provide appropriate expertise but it also introduces scien- 
tists outside of NASA to NASA problems. A small group of core RIACS staff provides continuity and inter- 
acts with an ARC technical monitor and scientific advisory group to determine the RIACS mission. RIACS 
activities are reviewed and monitored by a USRA advisory council and ARC technical monitor. 

Research at RIACS is currently being done in the following areas: 

Parallel Computing 

Advanced Methods for Scientific Computing 
Learning Systems 

High Performance Networks and Technology 
Graphics, Visualization and Virtual Environments 

In the past year, parallel compiler techniques and adaptive numerical methods for flows in complicated 
geometries were identified as important problems to investigate for ARC’S involvement in the Computa- 
tional Grand Challenges of the next decade. One of RIACS stafF scientists. Dr. Robert Schreiber, has 
devoted most of his time to parallel compilers, and a senior visitor. Professor Marsha Berger, who works 
on adaptive grid methods, was invited to visit RIACS from September 1, 1991 through August 31, 1992. 

During the past six months Professor Petter Bjorstad, of the University of Bergen, has begun a visit to 
RIACS (7/1/92 - present). He will work on domain decomposition techniques for PDEs for use on parallel 
computers. Dr. Pelle Olsson, from the University of Uppsala, joined RIACS as a post doctoral scientist in 
November 1993. He is working on the derivation of boundary conditions for the Navier-Stokes equations 
which are stable for nonsmooth boundaries and their efficient implementation on parallel computers. 

We concluded a summer student visitors program during this six months. We had six visiting graduate stu- 
dents that worked on projects over the summer and presented seminars on their work at the conclusion of 
their visits. 

RIACS technical reports are usually preprints of manuscripts that have been submitted to research journals 
or conference proceedings. A list of these reports for the period July 1, 1992 through December 31, 1992 
is in the Reports and Abstracts section of this report. 
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II. RESEARCH IN PROGRESS 


A. PARALLEL COMPUTING 


HIGH PERFORMANCE FORTRAN FORUM 

Robert S. Schreiber 

Schreiber is an active participant in the High Performance Fortran Forum (HPFF), an academic, govern- 
ment, industry group that is designing a portable language based on Fortran 90 for programming distributed 
memory parallel machines. The Forum has approximately 40 participants. The group has met every six 
weeks throughout 1992 in order to complete its work and present a draft of the language at the Supercom- 
puting ‘92 Conference. Schreiber has been working on the issues of language compatibility with Fortran 
77, on the semantics of the data distribution directives, and on the library of new intrinsic functions. He has 
been chairman of the Intrinsics Committee. 

HPFF is completing its work on schedule. A draft version was presented at Supercomputing ‘92 and is now 
available for public comment via ftp. A final meeting in March 1993 will be used to finish the final draft. 

A number of hardware and software vendors have announced their intention to develop and market HPF 
compilers. It appears that HPF will be accepted by the parallel computer industry and become a de facto 
standard. If this happens, the Forum’s work will have been a success. 


RESEARCH IN HIGHLY PARALLEL MATRIX COMPUTATION 
Robert S. Schreiber 

Robert S. Schreiber, with Sid Chatteijee and Edward Rothberg, is continuing to investigate scalable meth- 
ods for sparse Cholesky factorization on distributed memory machines. These methods must use block dis- 
tribution (rather than row or column distribution) of the sparse matrices. Load balance, communication 
load, and scheduling are all important issues in achieving high performance on highly parallel machines. 

With Elizabeth Jessup, they are also looking into scalable methods for the dense, symmetric eigenvalue 
problem. To date, there has been no satisfactory method implemented for the reduction of such a matrix to 
tridiagonal form with orthogonal similarity transformations. Such a method could be the basis for a stan- 
dard highly parallel eigenvalue solver, since parallel methods for the tridiagonal eigenvalue problem are 
now available. 

With Petter Bjorstad, a software message routing system for grid and toroidal machines has been devel- 
oped in order to assess their usability for solving PDEs with unstructured meshes. This work will continue 
during the coming six-month period. 

Also, in progress: Efficient, Massively Parallel Eigenvalue Computation, with Yan Huo (Electrical Engi- 
neering Department, Princeton University). 
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MASSIVELY PARALLEL COMPUTERS 
AND COMPUTATIONAL ELECTROMAGNETICS 
Niel K. Madsen 

Madsen has completed the parallelization of the preprocessor, PREDSI3D, for the time domain electromag- 
netic modeling code DSI3D. Initial performance tests have indicated that some parts of the preprocessor 
perform well on the Intel iPSC/860 computer, while other parts perform less well. A number of the subrou- 
tines in PREDSI3D are required to perform disk read and write operations and we have found that disk i/o 
can be a real bottleneck on this particular computer. As there are only ten i/o nodes, when more than ten 
nodes need to perform i/o linear speedup with the number of processors does not occur. Methods are cur- 
rently being investigated to minimize the i/o problems. 

The basic data allocation scheme used to distribute the work among the processors has been changed. Pre- 
viously, the actual problem field component variables were distributed. This led to the need to solve very 
large partitioning problems. Now the problem grid cells are distributed among the processors. As there are 
as many as six field component variables associated with each cell, this strategy reduces the magnitude of 
the partitioning problem by a corresponding factor of six. This change of partitioning strategy required rel- 
atively minor changes to the preprocessor. 

As replacements for the existing parallel machines are imminent, the DSI3D codes will be ported to these 
new machines - which will in principle provide a significantly greater 3D modeling capability. As in the 
past, it is expected to challenge these new machines with unstructured grid electromagnetic problems 
which are larger than have been solved heretofore. 


B. ADVANCED METHODS FOR SCIENTIFIC COMPUTING 

ITERATIVE METHODS FOR NON-HERMTTIAN LINEAR SYSTEMS 

Noel M. Nachtigal 

Nachtigal, together with Roland Freund at AT&T Bell Labs, continued work in the area of iterative methods 
for non-Hermitian linear systems. In particular, they were interested in developing an alternate implemen- 
tation of the quasi-minimal residual (QMR) method, based on a coupled two-term recurrence version of the 
underlying nonsymmetric Lanczos algorithm, rather than the classical three-term Lanczos recurrence. The 
motivation behind this work is that, in finite precision arithmetic, algorithms based on coupled two-term 
recurrences seem to have better numerical properties than the mathematically equivalent algorithms based 
on three-term recurrences. 

They succeeded in developing a coupled version of the Lanczos process, complete with a look-ahead strat- 
egy that is necessary to avoid breakdowns. They have then applied the quasi-minimal residual approach, 
obtaining a new implementation of the QMR algorithm. Preliminary numerical experiments indicate that 
the new version does indeed have numerical properties that are superior to those of the original QMR algo- 
rithm. 

Work is in progress on writing a software package that implements the new version of the QMR algorithm 
for the solution of non-Hermitian linear systems. When finished, it is planned to make this software avail- 
able to the scientific community, as has been done in the past with the original QMR implementation. As 
an aside, Thinking Machines Corporation has ported the codes made available with the original QMR algo- 
rithm to the Connection Machine and has incorporated this version in their CM Scientific Software Library. 
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and refinement algorithm is the key to the success of this procedure. The data structure for this algorithm 
is implemented in C and consists of a series of linked lists. These linked lists allow the mesh connectivity 
to be rapidly reconstructed when individual mesh points are added and/or deleted. It also allows for aniso- 
tropic refinement of the mesh. 

Aerodynamics calculations performed on structured-grids have difficulties resolving localized flowfield 
features such as shocks, vortices, and acoustic waves. Unstructured-grid models can make use of localized 
mesh refinement to resolve these flow features. However, this mesh refinement is only effective if it can be 
performed efficiently in three dimensions. This new procedure for dynamic mesh adaption directly 
addresses this problem by using an innovative data structure that is well suited for large-scale computa- 
tions. When coupled with a 3-D unstructured-grid Euler solver, the mesh adaption scheme will provide 
accurate solutions for complex aerodynamic flowfields. 

The dynamic mesh adaption scheme is currently being tested for large problems on a Cray Y-MP computer. 
Particular attention is focused on CPU time and memory requirements. Future work will test the perfor- 
mance of the mesh adaption scheme and 3-D Euler solver on a massively parallel computer system. 


RESEARCH IN DOMAIN DECOMPOSITION 

Petter Bjorstad 
Barry Smith (UCLA) 

This work consists of two main lines of development. Petter Bjorstad and Barry Smith (UCLA), are writing 
up a more complete description and discussion of various algorithms that have been published over the 
years (this work takes the form of a book manuscript). A significant part of this work is to recast algo- 
rithms into a consistent form and describe them such that their implementation on parallel computers 
becomes feasible for non-experts in the field. This is important in order to get this powerful class of algo- 
rithms used in real application work. Along with this work is also actual software engineering and parallel 
im plementations, in particular targeted towards problems from CFD calculations (this research is also joint 
with Eric Grosse and Bill Coughran of Bell Labs). 

State of the art iterative methods for the solution of large scale discrete systems of equations arising from 
the discretization of partial differential equations are the so called Krylov Space Methods. These methods 
require the efficient application of a preconditioner as well as the discrete operator acting on a vector. 
Large scale computations in 3 dimensions often make use of highly unstructured grids. It is, therefore, of 
interest to investigate if the necessary communication in a parallel computer can be carried out without a 
large extra cost in this case. This research is joint with Robert Schreiber and considers the development of 
a communication compiler for 2-D, mesh interconnected parallel computers. In particular, data parallel 
(SIMD-style) machines like the MasPar, are considered. 

RESEARCH IN HIGHLY PARALLEL COMPUTATION 
Petter Bjorstad 
Rob Schreiber 
Erik Boman (Stanford) 

Petter Bjorstad and Rob Schreiber, have developed a software message routing system for grid and toroidal 
machines, in order to assess their usability for solving PDEs with unstructured meshes. They will be con- 
tinuing this work during the coming six-month period. 

Bjorstad, with Erik Boman and Rob Schreiber, has also initiated the development of a high performance, 
modular implementation of the important BLAS-module DGEMM for parallel, distributed memory comput- 
ers. 
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images of Ganymede, where the curvature of the planetary surface must be taken into account. 

The development of a procedure for combining multiple image information required the solution to a num- 
ber of difficult problems. It was particularly important to find very accurate registrations of the individual 
images to the composite, correct for nonuniform gain in the camera, learn and correct for the point-speed 
function within the camera, and other corrections related to noisy pixels. Algorithms were found for solv- 
ing these problems and also for taking into account the neighbor correlations. The results of this work are 
reported in a NASA tech note, and have been submitted for presentation at a conference. 


AUTOCLASS - AN AUTOMATIC CLASSIFICATION SYSTEM 
Peter Cheeseman 
John Stutz (NASA Code FIA) 

Robin Hanson (NASA/Sterling) 

Will Taylor (NASA/Sterling) 

For several years, this group has been developing probabilistic techniques for the automatic classification 
of data. This is a useful technique for exploratory data analysis. This research has been implemented and 
tested in the AutoClass program. AutoClass accepts as input a set of independent data cases, each described 
by the same set of attributes. Then, using appropriate prior expectations and a likelihood model describing 
possible attribute relationships, AutoClass seeks the most probable classification given the data. It does not 
need to be told how many classes are present-it infers this from the data. 

Recent AutoClass developments include models for full and partial correlation, hierarchical interclass rela- 
tionships, and angular data. These developments, and the mathematics behind them, are described in a 
detailed NASA Tech Report and a conference paper. Copies of this experimental AutoClass version have 
been distributed to a number of researchers. A production version is currently distributed by NASA’s COS- 
MIC software center. A massively parallel version is now running on the NAS Connection Machine. Auto- 
Class continues to provide a useful tool both for internal use, and externally. 

AutoClass was applied to LandSat image data. It was found that AutoClass gives reasonable looking clas- 
sifications and that some classes have obvious interpretations. A NASA Tech Report was prepared describ- 
ing AutoClass and a classification of an intensively studied area in Kansas (First ISILCP Field Experiment), 
and collaborators familiar with the area are being sought. This LandSat work did not use spatial correla- 
tions among adjacent pixels and was unable to integrate information from successive images. These limi- 
tations led to the most recent project-surface reconstruction from multiple images. 

SPARSE DISTRIBUTED MEMORY 
Pentti Kanerva 

This is long-term basic research into the engineering principles by which the brain and the rest of the ner- 
vous system organize and process information. Uncovering these principles serves dual purposes: Build- 
ing of machines with skills, and designing of systems to fit human operators. The resulting technology will 
have wide use in aviation and in manned and unmanned space exploration. 

Sparse distributed memory provides natural entry into the research. It was developed as a mathematical 
and engineering model of human long-term memory, and it is at once an associative memory, a random- 
access memory, a sequence memory, an artificial neural network, and a model of the cerebellum. The 
memory stores long (e.g., 1,000-bit) words that need not be precise, and it is also addressed by such words. 
An object or a scene or a moment of experience can be represented in the memory by a single word that 
contains a large amount of information about the object or the scene or the moment, and the word can serve 
equally as data and as a memory pointer. Memory pointers with high information content are at least partly 
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high data-rate requirements and also because the project office is conveniently located at NASA Ames. 

During the past six months M. Johnson has developed a model of the end-to-end Space Station communi- 
cations system, using the simulation package Network II.5. The performance of the end-to-end system is 
dependent on such factors as network protocols, scheduling restrictions, limitations on uplink bandwidth, 
commanding procedure, space-to-ground communication outages, anticipated level of interaction between 
a scientist and his spacebome experiment, and safety considerations. During the next six months, M. 
Johnson will refine the model and use it to determine the effects of the above factors on system perfor- 
mance. 

The sponsor for this activity is Space Station Headquarters. The work is being conducted in collaboration 
with the Spacecraft Data Systems Research Branch of the Information Sciences Division, which is 
involved in analyzing broader data-management issues and conducting testbed activities in support of 
Space Station Freedom. Information about the CFP is being provided by project personnel. 

RIACS NATIONAL SOFTWARE EXCHANGE (NSE) PROJECT 

Mike Raugh 

The RIACS NSE Project continues supporting the HPCC Software Exchange. The Software Exchange (SE) 
Project, led by Barry Jacobs at Goddard Space Flight Center is near the end of its second year. Participants 
representing the major HPCC Program agencies - DARPA, DOC/NIST, DOE, and NSF, as well as NASA - are 
cooperating to make information of interest and relevance to HPCC Program software developers available 
on-line through a shared SE infrastructure. RIACS is working with Jacobs to further develop and apply the 
GSFC Logical Library and Logical Books systems as prototypes for an on-line library and an experimental 
uniform interface, respectively. The Logical Library and its main catalogue permit easy lookup and access 
of major information resources on the Internet, such as RECON, Netlib, GAMS, WAIS, Cugdus, and CASRD, 
to name a few. Approximately forty-five such resources are available at this time, about a dozen of which 
are available in the “books” format. 

The SE Project, and the RIACS supporting role, undergoes continual revision as strategies evolve rapidly 
for improving and populating the Library. The basic structure of the Library and Books systems have 
begun to stabilize, while the tools themselves are developing at a quickening pace. Beta versions of client/ 
servers for both Books and the Library are nearing completion and multimedia capabilities have been 
incorporated and demonstrated. The assigned role for RIACS in this period has been to continue work on 
Books as a uniform access mechanism, in particular to develop the initial version of the Library main cata- 
logue and publish it as a Book, to continue developing tools ingesting published materials into the Books 
format, and to develop tools for Library maintenance functions. Once these tools themselves have stabi- 
lized, they will be reported and made available to other HPCCP scientists and programmers who wish to 
publish Books. 

The primary objective for RIACS in this half-year has been to work with the DAVID Publishing Group at 
GSFC to refine the process of publishing on-line Books, to seek additional resources that are relevant to 
workers in the HPCC, and to develop a first edition of the Library’s catalogue of major resources. The term 
“Books” in this context refers to a particular way of representing information in both VT100 and, most 
impressively, X-Windows formats. 

Such Books are proving to be a versatile and powerful means of representing widely differing categories of 
data in a consistent format, including software, textual information, images, recorded sound, and bibliogra- 
phies. Books are being used even to describe hierarchies of information, leading from the more general 
levels of a taxonomy to more particular levels. A secondary objective for RIACS throughout the project has 
been to develop software for Books publishing. 
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1. DISTRIBUTED PROCESSING: . __ 

a) A master process runs on a host workstation running an XI 1 server. This process initiates the slave and 
provides a medium speed operator’s user interface to the slave process via a socket connection. 

b) A slave process runs on the high performance SGI SkyWriter that polls time critical user input devices 
(head tracker and 6 DOF joystick), responds to operator messages sent via the socket connection, and con- 
tinuously draws animation frames. This process provides the rapid (10 to 30 Hz) video update for viewing 
in a stereoscopic, head-mounted display. 


2. CONCURRENT RENDER PROCESSING: 

The SkyWriter rendering operation uses 3 concurrent processes (application, cull, and draw) running on up 
to 4 CPU’s. Each animation frame is pipelined through the three stages such that frame updates (the draw 
process) always occur as fast as possible with no waiting between updates. In dual channel stereo mode, 
the cull, draw, graphics hardware, and video output pipeline for each channel run in parallel. 

The VPE software was extended to provide multi-resolution terrain modeling using several level-of-detail 
(LOD) modeling criteria. For each region of terrain several polygon mesh models at many different sam- 
pling resolutions are stored in a quad tree. At each update frame analysis of each subregion of the terrain 
surface determines an appropriate model resolution to use when drawing the terrain. The LOD analysis 
exploits the quad tree hierarchy and employs three criteria for selecting the appropriate model: range from 
the viewpoint, angular distance from the center of the field-of-view, and a user specified metric for applica- 
tion related regions-of-interest. The LOD analysis considers all three factors weighted according to user 
specifications. The analysis/feedback loop can be adjusted to respond to system load and user preferences. 
This technique provides both a high frame update rate and high level-of-detail for terrain regions important 
to the user (close to the viewpoint, near the center of field of view, and areas critical to the application) at 
the sacrifice of low detail in regions of low interest for visual perception and application requirements. 
The increased performance using this technique has greatly enhanced the VPE system’s ability to provide 
compelling virtual environment simulations. This result of this work is being presented in a conference 
proceedings in February 1993. 

This project is sponsored by the Human Interface Research Branch (Code FLM) of the Aerospace Human 
Factors Division (Code FL), NASA Ames. Dr. Michael W. McGreevy, Code FLM, is the NASA collaborator. 


F. OTHER PROJECTS 

GLOBAL CHANGE RESEARCH 
Richard Johnson 

The Earth Sciences Division (ESD) at NASA Ames Research Center has been exploring the research options 
that could be appropriate for support by the NASA Office of Commercial Programs (OCP) and various pri- 
vate sector industries. Support has been provided to ESD for the development of an overall strategic plan 
and a set of project options for discussions with OCP. This plan and project options were presented to and 
discussed with OCP program managers at NASA Headquarters via a video conference in July. Two follow- 
up presentations to OCP personnel were made in August at NASA Ames. R. Johnson discussed the opportu- 
nities for research and commercialization in the Biogeochemical Dynamics research and technology arena. 
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III. TECHNICAL REPORTS 


AMR on the CM-2 

MARSHA J. BERGER 

JEFF SALTZMAN (Los Alamos National Laboratory) 

TR 92.16 August 1992 (17 pages) 

Submitted to Applied Numerical Mathematics 

We describe the development of a structured adaptive mesh algorithm (AMR) for the CM-2. We develop a 
data layout scheme that preserves coarse grids. On 8K of a 32K machine we achieve performance slightly 
less than 1 CPU of the Cray Y-MP. We apply our algorithm to an inviscid compressible flow problem. 

Optimal Evaluation of Array Expressions on Massively Parallel Machines 
SIDDHARTHA CHATTERJEE 

JOHN R. GILBERT (Xerox Palo Alto Research Center) 

ROBERT S. SCHREIBER 

SHANG-HUA TENG (Xerox Palo Alto Research Center) 

TR 92.17 September 1992 (36 pages) 

To appear in the Proceedings of the Second Workshop on Languages, Compilers and Runtime Environments for Dis- 
tributed Memory Multiprocessors, Boulder, CO., September 30 - October 2, 1992. 

We investigate the problem of evaluating Fortran 90-style array expressions on a massively parallel distrib- 
uted-memory machine. On such machines, an elementwise operation can be performed in unit time for 
arrays whose corresponding elements are in the same processor. If the arrays are not aligned in this man- 
ner, the cost of aligning them is part of the cost of evaluating the expression. The choice of where to per- 
form the operation then affects this cost. We present a dynamic programming technique to solve this 
problem efficiently for a wide variety of interconnection schemes, including multidimensional grids and 
rings, hypercubes and fat-trees. We also consider expressions containing operations that change the shape 
of the arrays, and show that our approach extends naturally to handle this case. 


Automatic Array Alignment in Data-Parallel Programs 
SIDDHARTHA CHATTERJEE 

JOHN R. GILBERT (Xerox Palo Alto Research Center) 

ROBERT S. SCHREIBER 

SHANG-HUA TENG (Xerox Palo Alto Research Center) 

TR 92.18 October 1992 (13 pages) 

To appear in the Proceedings of the Twentieth Annual ACM SIGACT/SIGPLAN Symposium on Principles of Program- 
ming Languages, Charleston, SC, January 1993. 

Dataparallel languages like Fortran 90 express parallelism in the form of operations on data aggregates 
such as arrays. Misalignment of the operands of an array operation can reduce program performance on a 
distributed-memory parallel machine by requiring nonlocal data accesses. Determining array alignments 
that reduce communication is therefore a key issue in compiling such languages. 

We present a framework for the automatic determination of array alignments in dataparallel languages 
such as Fortran 90. Our language model handles array sectioning, reductions, spreads, transpositions, and 
asked operations. We decompose alignment functions into three constituents, axis, stride, and offset. For 
each of these subproblems, we show how to solve the alignment problem for a basic block of code, possi- 
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IV. PUBLICATIONS 


RUPAK BISWAS 

Applications of an Unstructured Adaptive-Grid Euler Solver to Helicopter Aerodynamics and Acoustics, 
R. C. Strawn, R. Biswas, and M. Garceau, Compendium of Abstracts of the NASA Computational Aero- 
sciences (CAS) Conference, NASA Ames Research Center, Moffett Field, CA, August 18-20, 1992, pp. 66-67. 

A New Procedure for Dynamic Adaption of Three-Dimensional Unstructured Grids, R. Biswas and R. C. 
Strawn, to appear in the Proceedings of the 31st AIAA Aerospace Sciences Meeting, Reno, NV, January 11- 
14, 1993. 

PETTER BJORSTAD 

Two Different Data-Parallel Implementations of the BLAS (with Tor Sorevik, University of Bergen). To 
appear in NATO proceedings from a meeting on parallel computing, Cetraro, Italy, June. 

Dataparallel BLAS As a Basis for LAPACK on Massively Parallel Computers (with Tor Sorevik, University 
of Bergen). In NATO Advanced Study Institute , published by Kleuver 1993. 

Large Scale Structural Analysis on Massively Parallel Computers (with Jeremy Cook, University of Ber- 
gen). In NATO Advanced Study Institute , published by Kleuver 1993. 

SIDDHARTHA CHATTERJEE 

Optimal Evaluation of Array Expressions on Massively Parallel Machines, by Siddhartha Chatteijee, John 
R. Gilbert (Xerox PARC), Robert Schreiber, and Shang-Hua Teng (Xerox PARC). Appears in the Proceed- 
ings of the Second Workshop on Languages, Compilers, and Runtime Environments for Distributed Mem- 
ory Multiprocessors, Boulder, CO, September 1992. Also available as RIACS TR 92.17. 

Automatic Array Alignment in Data-Parallel Programs, by Siddhartha Chatterjee, John R. Gilbert (Xerox 
PARC), Robert Schreiber, and Shang-Hua Teng (Xerox PARC). To appear in the Proceedings of the Twenti- 
eth Annual SIGPLAN-SIGACT Symposium on Principles of Programming Languages, Charleston, SC, Janu- 
ary 1993. Also available as RIACS TR 92.18. 

Compiling Nested Data-Parallel Programs for Shared-Memory Multiprocessors, by Siddhartha Chatteijee. 
Accepted for publication in ACM Transactions on Programming Languages and Systems. 

Generating Local Addresses and Communication Sets for Data-Parallel Programs, by Siddhartha Chatter- 
jee, John R. Gilbert (Xerox PARC), Fred J. E. Long (UC Santa Cruz), Robert Schreiber, and Shang-Hua Teng 
(Xerox PARC). To appear in the Proceedings of the Fourth ACM SIGPLAN Symposium on Principles and 
Practice of Parallel Programming, San Diego, CA. 

Implementation of a Portable Nested Data-Parallel Language, by Guy E. Blelloch (CMU), Siddhartha Chat- 
teijee, Jonathan Hardwick (CMU), Jay Sipelstein (CMU), and Marco Zagha (CMU). Tb appear in the Pro- 
ceedings of the Fourth ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming, San 
Diego, CA. 

LEWIS E. HITCHNER 

Virtual Planetary Exploration: A Very Large Virtual Environment, ACM SIGGRAPH ’92, course notes for 
tutorial on “Implementing Immersive Virtual Environments,” Steve Brsyon, course chair, Chicago, IL; 
August 1992. 

The NASA Ames Virtual Planetary Exploration Testbed, Proceedings of IEEE Wescon ’92, Anaheim, CA, 
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V. RIACS SEMINARS AND COLLOQUIA 


Justin D. Paola - RIACS/NASA Ames Research Center - University of Arizona 
The Application of Neural Networks in the Classification of Multi-Spectral 
Satellite Imagery 

Shang-Hua Teng - XEROX Palo Alto Research Center 
A Geometric Approach to Mesh Partitioning 

Marsha Berger - RIACS/NASA Ames Research Center - New York University 
Issues for Adaptive Mesh Refinement: Parallelism and Geometry 

Petter Bjorstad - RIACS/NASA Ames Research Center -University of Bergen 

Benchmarking of Scientific Computers, SLALOM - A Case Study 

James M. Strzelec - RIACS/NASA Ames - Stanford University 
Iterative Methods: An Informal Comparison 

Bruce Hendrickson - Sandia National Laboratories 
Spectral Octasectionfor Mapping Parallel Computations 

Mike Raugh - RIACS/NASA Ames Research Center 

Towards an HPCC Internet Library - An Informal Discussion of the 

HPCC Software Exchange Prototype 




July 29 


August 1 1 


August 25 


Sept. 8 


Sept. 9 


October 26 


October 30 
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LEWIS E. HITCHNER 

Co-taught a tutorial on Implementing Immersive Virtual Environments at the annual ACM SIGGRAPH con- 
ference, Chicago, IL, August, 1992. 

Visited the Alpha_l Graphics Research Group, Computer Science Dept., Univ. of Utah, and presented a sem- 
inar on the NASA Ames virtual environment research, Salt Lake City, Utah, September 1992. 

MARJORY J. JOHNSON 

Member of the program committee and session chair for the 4th IFIP Conference on High Performance Net- 
working, Liege, Belgium, December 1992. 

PENTTI KANERVA 

Plenary talk on Associative-Memory Models of the Cerebellum at the Second International Conference on 
Artificial Neural Networks (ICANN’92) in Brighton, UK, September 3-7. 

NOEL M. NACHTIGAL 

Presented An Implementation of the QMR Method Based on Coupled Two-Term Recurrences, for a confer- 
ence at the NATO Advanced Study Institute, Leuven, Belgium, August 3-14, 1992. 

Presented An Implementation of the QMR Method Based on Coupled Two-Term Recurrences, at the SIAM 
40th Anniversary Meeting, Los Angeles, CA, July 20-24, 1992. 

Presented Iterative Methods in Linear Algebra, at the Naval Postgraduate School, Monterey, CA, November 
19, 1992. 

Presented An Implementation of the QMR Method Based on Coupled Two-Term Recurrences, at the Uni- 
versity of Karlsruhe, Germany, September 11, 1992. 

Presented a seminar, Iterative Methods for Non-Hermitian Linear Systems, at the IBM Scientific Center, 
Heidelberg, Germany, September 10, 1992. 

Participated in the five-day course on Large Scale Scientific Computation at the University of Bielefeld, 
Bielefeld, Germany, August 31 - September 4, 1992. 

JOSEPH OLIGER 

Elected to the Society of Industrial and Applied Mathematics Board of Trustees, 1993-1996. 

Chairman, National Center for Atmospheric Research Scientific Computing Division Advisory Commit- 
tee. 

Member, Organizing Committee for IMA Summer Program on Modeling, Mesh Generation and Adaptive 
Numerical Methods for Partial Differential Equations, 5-23 July 1993. 

ROBERT S. SCHREEBER 

Presented High Performance Fortran Intrinsics at Workshop on High Performance Fortran, Supercomput- 
ing ‘92, November 17, 1992. 

Presented Are Sparse Matrices Poisonous for Highly Parallel Machines for UNI-C, (Danish Supercomputer 
Center) Copenhagen, Denmark, July 1992. 

Presented Language and Compiler Grand Challenges, an invited talk; Gordon Research Conference on 
Software Tools and Libraries for Concurrent Supercomputing, Plymouth, NH, July 27-3 1 , 1992. 
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VII. RIACS STAFF 


ADMINISTRAT IVE STAFF 

Joseph Oliger, Director - Ph.D., Computer Science, University of Uppsala, Sweden, 1973. Numerical 
Methods for Partial Differential Equations (03/25/91 - present). 

Frances B. Abel, Office and Financial Manager (5/5/88 - present). 

Deanna M. Gearhart, Administrative Assistant II (5/9/88 - present). 

Anne Kohutanycz, Administrative and Systems Assistant (1/21/85 - 8/20/92). 

Evangeline Tanner, Administrative Assistant I (4/5/90 - present). 

TECHNICAL S UPPORT STAFF 

Bradley Christofferson, A.A. - Computer Information Systems, DeAnza College 1987. Computer Sys- 
tems Manager (10/23/89 - 10/30/92) 

Jessica Casillas, B.A. - Mathematics, San Francisco State University 1986. Computer Systems Manager 
(11/25/92 - present). 

RIACS SCIENCE COUNCIL 

Dr. David Cummings (Interim Convener), Executive Director, Universities Space Research Associa- 
tion, Columbia, MD. 

Dr. Dennis B. Gannon, (Covener), Director, Center for Innovative Computer Applications (CICA), Indi- 
ana University, Bloomington, IN. 

Dr. Joseph Flaherty, Chair, Department of Computer Science, Rensselaer Polytechnic Institute, Troy, NY. 

Dr. Joseph Oliger (Ex-Officio), Director, Research Institute for Advanced Computer Science, NASA 
Ames Research Center, Moffett Field, CA. 

Dr. James W. Demmel, Computer Science Division, University of California, Berkeley, CA. 

Dr. David Gottlieb, Division of Applied Mathematics, Brown University, Providence, RI. 

Dr. Kenneth W. Neves, Boeing Company, Seattle, WA. 

Dr. Thomas H. Pulliam, NASA Ames Research Center, Moffett Field, CA. 

Dr. Daniel A. Reed, Department of Computer Science, University of Illinois, Urbana, IL. 

Dr. Robert B. Schnabel, Department of Computer Science, University of Colorado, Boulder, CO. 

Dr. Marc Snir, IBM Thomas J. Watson Research Center, Yorktown Heights, NY. 
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parabolic PDEs, numerical methods for PDEs on parallel computers (11/2/92 - present). 

STUDENTS 

Fred Long - Computer Science, UC San Diego, (1990). Construction, implementation and evaluation of 
algorithms for efficient address generation for FORTRAN 90 programs mapped onto distributed-memory 
parallel computers (6/22/92 - 8/28/92). 

Justin Paola - Electrical Engineering & Computer Science, UC Berkeley, (1990). Image processing/ 
remote sensing (5/26/92 - 8/7/92). 

Brian Rogoff- Mathematics, M.I.T. (1986). Scientific visualization (7/1/92 - 9/31/92). 

James M. Strzelec - Applied Mathematics, SUNY, Stony Brook, (1992). Iterative methods (6/22/92 - 
9/9/92). 

Steven Suhr - Computer Science, Stanford University, (1987). Algorithms for scientific computing; pro- 
gramming languages (7/1/92 - same). 

CONSULTANTS 

Marsha Berger - New York University. Computational fluid dynamics; parallel computing (9/1/92 - 
present). 

Tony F. Chan - Professor of Mathematics, University of California, Los Angeles. Efficient algorithms in 
large-scale scientific computing, parallel algorithms and computational fluid dynamics (10/01/86 - 
present). 

Michael Eiermann - University Karlsruhe, Germany. Iterative methods in Linear Algebra (9/1/92 - 
12/31/92). 

Paul O. Frederickson - CRAY Research. Parallel algorithms, with emphasis on multigrid algorithms and 
their application to fluid dynamics and particle models (9/25/91 - 12/31/92). 

John Gilbert - Research Scientist, Xerox Palo Alto Research Center. Parallel computing and theoretical 
computer science (5/1/92 - present). 

Elizabeth Jessup - Professor, University of Colorado, Boulder. Parallel matrix computation; eigenvalue 
problems; programming distributed memory systems (10/19/92 - 12/31/92). 

Niel K. Madsen - Lawrence Livermore National Laboratory. Numerical solutions of partial differential 
equations, with specific interests in method of lines techniques, PDE software, matrix algorithms for vector 
and parallel computers (10/18/90 - present). 

Jeffrey D. McDonald - Applications Specialist, MASPAR Computer Corporation. Computational fluid 
dynamics (5/1/92 - present). 

Alex Pothen - Professor of Mathematics, University of Waterloo, Ontario, Canada. Parallel scalable 
sparse matrix algorithms (10/01/92 - present). 

David Rogers - Project Manager, Molecular Simulations, Inc. Next-generation computational chemistry 
software product (2/28/92 - 12/31/92). 
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